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Abstract 

We study phenomenological aspects of supersymmetric SO(W) 
GUTs with sum rules among soft SUSY breaking parameters. In 
particular, the sum rule related to the stau mass leads to the con- 
straints from the requirements of successful electroweak breaking and 
the positivity of stau mass squared. The bottom quark mass is also 
estimated. 
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The unification of force based on the minimal supersymmetric standard 
model (MSSM) |lj has been hopeful from the data of precision measurements 
[0]. The supersymmetric SO (10) grand unified theory (SUSY SO (10) GUT) 

is one of attractive candidates of realistic theory at and above the GUT 
scale Mx because a simplest unification of quarks and leptons can be realized 
in each family. After the gauge symmetry breakdown, the remnant exists in 
the MSSM as specific relations among physical parameters at Mx, which are 
usually used as initial conditions on the analysis by the use of renormalization 
group equations (RGEs), e.g., g$ = g% = g\ = g for the gauge couplings of the 
SM gauge group G SM = SU(3) C x SU(2) L x U(1) Y and = Yf = Yf = 
yafi f Qr Yukawa coupling matrices of up-type quarks, down-type quarks 
and leptons. Here a, j3 are the family indices. The number of independent 
soft SUSY breaking parameters is also reduced by the SO (10) symmetry, i.e., 
the parameters are given as (m^, tuh, M, A a @, B) at Mx for soft squark 
and slepton masses, soft Higgs masses, gaugino masses, A-parameters and 
Z?-parameter up to the contributions of SO (10) breaking including D-term 
contribution ||, ||. 

The magnitude of each parameter is expected to be determined by some 
underlying theory or new concept of SUSY-GUT. The supergravity (SUGRA) 
is regarded as a promising theory, describes the physics beyond SUSY-GUT 
effectively, and offers an interesting scenario of the origin of soft SUSY break- 
ing parameters 0. The structure of the SUGRA model is reflected on the 
pattern of the parameters, e.g., the universal soft SUSY breaking parameters 
originate in models with a canonical Kahler potential. The analyses based 
on this type of initial conditions have been intensively carried out [[/J. 

There is another interesting scenario to control the relations among phys- 
ical parameters. The parameters can be reduced by the adoption of a new 
concept 'coupling reduction' by the use of RG invariant relations [§]. The 
assumption in a model is that the Yukawa couplings are expressed RG- 
invariantly by the gauge coupling g as 

Yijk — QPijkj (1) 

where are model-dependent constant independent of g at the tree level. 
Higher order corrections are systematically calculated. It is called gauge- 
Yukawa unified (GYU) model 0. Recently, it is found that the following 
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relations of soft parameters are RG-invariant [TO , |TT|l,fl 



m 2 = m 2 + m) + m\ = M 2 , (2) 
= —M (3) 

where 3^*., rrii and are Yukawa couplings, soft scalar masses and A- 
parameters. The indices k denote the particle species. It has been 
known that the relations (0) and (§) are derived in other several theories, 
i.e., the finite field theories |]I3"|| , a certain class of 4-dimensional models from 
superstring theory (SST) fl4| and a non-minimal SUGRA model with a cer- 
tain type of structure regarded as a generalization of models from SST [[llj . 
Hence this type of relations can give a hint to high energy physics beyond the 
MSSM. Therefore it is important to study the phenomenological implications 
and low-energy consequences from this type of relations [15|, [16| . 

The low-energy constraints from (0) and ([]) are studied in Ref. []1]J based 



on finite SUSY SU(5) GUT models and it is shown that eigenvalues of stau 
masses squared m~ tend to be negative in some parameter region. The 
sbottom and stop fields always have heavier masses than the lightest stau. 
Thus the condition of the positivity m~ > as well as the electroweak 
symmetry breaking conditions constrains the parameter regions severely in 
this type of models. 

In the sector where the above relations (|T|), (Q) and (|3|) hold on, indepen- 
dent parameters are limited to g, M and mj. The radiative corrections of 
m 2 are given as functions of g and M because the contribution from Yukawa 
couplings contains soft scalar masses only as a combination of YJm 2 . In this 
sense, this type of soft SUSY breaking parameters are much more restrictive 
than a general non-universal one. 

Let us compare the excluded regions of soft stau mass mf and gaugino 
mass M at Mx between the models with the universal type of soft masses 
(M, mo) inspired by a minimal SUGRA and the models with the relations (@) 
and (0) at Mx, based on the MSSM with large tan (3. From the requirement of 
successful electroweak symmetry breakdown, there is the constraint m 2 Hi — 
m 2 H2 > Mf at the weak scale. Here m 2 Hi and m 2 H are soft Higgs masses 
squared with the hypercharge —1/2 and 1/2, respectively. Using the analysis 
of RGEs, we get the relation such that m 2 Hi —m 2 H2 = C\M 2 — c 2 fn 2 /3 under 
the condition that rr? H = m 2 H2 at Mx [0 and the constant factors C\ and c 2 



We can calculate higher order corrections to RG invariant relations E 
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are of 0(0.2). From these formulae, the region with m?(= mo) > M, (m? ^> 
Mz) is excluded in the universal case.Q On the other hand, in the case with 
sum rules, the excluded region is M < O(200)GeV independent of the value 
of m f . Another requirement is the positivity of physical stau mass squared. 
The Yukawa coupling induces to a radiative correction with a negative sign 
to the mass squared. Thus the stau mass squared can be negative if the 
magnitude of YJ m2 is sizable in the large tan /3 scenario. This happens easily 
in the case with sum rules because of the existence of the relation YJ m2 — M 2 
and, in this case, the region with M ^> nif is excluded. In the universal 
case, the situation is different because the contribution including the factor 
Em 2 (= 3ml) becomes tiny when the value of mo = mf is small, i.e. M 3> 
m . Hence it is expected that the excluded regions of (M, m,f) are located at 
the opposite corners (besides M < O(200)GeV in the latter case) each other 
from the above two phenomenological requirements. 

In this paper, we study generic SUSY 5*0(10) GUT model with sum rules 
among soft SUSY breaking parameters and make sure the above estimation 
of the parameter regions, quantitatively, imposing the conditions of success- 
ful electroweak symmetry breaking and the positivity of physical stau mass 
squared. The method of analysis is almost same as that made in Ref. [II . 
SUSY corrections to the bottom quark mass are also estimated. 

First we give a brief review on the derivation of the relations ([I]), (0) 
and (H) in GYU-models. We assume that parameters Y^, mf and Ajfc are 
expressed in terms of g and M. The relations (0), (||) and (ff) are obtained 
by solving so-called reduction equations perturbatively, 



o a dm f j_ a dm f j_ a 9m i 



dAjjk dA ijk dA ijk 



where fix denotes a /3-function of parameter X. For the application on an 
explicit model, see Ref. [BO]. 



t For the study on case with non-universal initial conditions for soft SUSY breaking 



parameters, see Refs. |18, 191. 
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Second we give some basic assumptions and relations on our analysis. 
The first assumption is that all of quarks and leptons in each family belong 
to one 16-plet under SO (10) and this 16-plet has the Yukawa coupling such 
as (16) 2 H where H is 10-plet including Hi and H 2 . We ignore the family 
mixing effects. The second one is that the relations (0) and @ hold on at 
Mx in the third family. The third one is that there are no extra contributions 
on the symmetry breaking SO (10) — > Gsm- (Later we relax this assumption 
by the introduction of D-term contribution.) Our initial conditions at Mx 
are summarized as follows, 

9s = 92 = 9i = 9, Y t = Y b = Y T = Y, 
m i = m l =m t = = m k = m ?6> 

m2 H ± = m H 2 = m H, 

A t = A b = A T = —M, 2m 2 6 + m 2 H = M 2 (7) 

where Q, t, b, f L and f R denote the SU(2) L doublet squark of the third 
family, the singlet stop, the singlet sbottom, the stau in the doublet slepton 
of the third family and the singlet stau. Here and hereafter we omit the 
index representing the third family. 

Next we parametrize the Yukawa coupling Y using g as Y = pg. The value 
of p gives an important information on the matter content and its interactions 
in GYU-models and/or the structure of superpotential in SUGRA. We take 
the following input parameters, 

M T = 1.777GeV, M z = 91.188GeV, (8) 
a^(Mz) = 127.9 + ^ log ^, (9) 

97T Mz 

sin 2 6 W (M Z ) = 0.2319 - 3.03 x 10~ 5 T - 8.4 x lO^T 2 (10) 

where T = M t /[GeV] — 165. Here M T and M t are physical tau lepton and 
top quark masses. The Yukawa unification condition gives the predicted top 
quark mass from the above experimental value of M T for each value of p. Fig. 
1 shows the predicted value of the physical top quark mass for k = p 2 . Thus 
we find the realistic region such that 0.7 < k < 1.4 to obtain the present 
experimental value of the top mass, M t = 175.6 ± 5.5GeV. For example, the 
value k = 1.0 leads to M t = 175 GeV and tan/? = 53, while tan/? = 50 and 
55 for k — 0.7 and 1.4, respectively. 
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Fig. 1: The top mass M t for k. 

Similarly we can calculate the bottom quark mass at the tree level for 
each value of k. However, SUSY corrections to the bottom quark mass is 
sizable in the large tan (3 scenario [[21] and that leads to another constraint 
H], IT3]. Thus, we will estimate the predicted bottom quark mass with SUSY 



corrections after calculations of the SUSY mass spectrum. 

We determine the values of fi and S-parameters by using the following 
two minimization conditions of the Higgs potential at the weak scale, 

2jiB 



777.1 + 777 9 



771, 



sin 2/3' 
-cos2/3(M; 



+ 777^ + m^) 



(11) 

(12) 



777 



Hi,-f7 2 



+ at 



where m\ 2 

Soft SUSY breaking parameters can be constrained from requirements. 
One of the most important constraints is the realization of electroweak sym- 
metry breaking. To this end, the following condition should be satisfied, 



(13) 



In addition, the bounded-from-below condition along the D-flat direction in 
the Higgs potential requires 



m\ + m\ > 2\fiB\. 



(14) 
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Another important condition is the positivity of physical scalar mass squared 
0. For example, two stau masses squared, m| x and m? 2 are obtained as 
eigenvalues of the following (mass) 2 matrix: 

m l L + M|cos2/5(-i + sm 2 9 w ) vY T (A T cos (3 — fj, sin (3) \ . . 
wr T (A T cos/3-/isin/5) m\ R - M| cos 2/3 sin 2 W J ^ ' 

where v 2 = (H 2 ) 2 + (Hi) 2 . Here we neglect the SUSY stau mass squared M 2 . 

Actually, in Ref. [|ll|], it is shown these conditions constrain severely the 
parameter space in SU(5) models with large tan/3. Because, in the large 
tan/3 scenario, the stau mass squared receives as sizable negative corrections 
due to the large tau Yukawa coupling as the soft Higgs masses squared m 2 Hi 
and rn 2 H2 do. Large values of m~ and m~ at Mx are favorable to avoid 
m| 12 < 0. Here m fl denotes the lightest mass of them. For example, it is 
impossible to satisfy these conditions in explicit SU(5) models in a small 
value of M. It is shown that the case with a common soft scalar mass, 
m 2 = M 2 /3, is not allowed in some finite SU(5) models. 

Now we discuss these constraints in generic SO( 10) model. Fig. 2 shows 
excluded regions by these constraints for k = 1.0 (tan/3 = 53). In this 
figure, the dotted region in the left side denotes the region forbidden by the 
electroweak breaking conditions. On the other hand, the place with asterisks 
correspond to the region with m~ < and squares denote the region where 
the light stau mass squared m~ is smaller than the lightest neutralino mass 
squared m 2 . Note that, in the case with the initial condition m 2 = m 2 , 

the lightest stau f\ almost originates in tr. Because the mass squared m~ 
has sizable positive radiative corrections due to SU(2) gauginos, and a half 
size of negative contribution from r Yukawa coupling compared with that to 
m| fl . In the whole parameter space, sbottom and stop are heavier than the 
lightest stau. 

Ht would be necessary to calculate the decay rate to the unboundcd-from-below direc- 
tion, e.g. corresponding to < in order to exclude completely such parameter region 
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Fig.2: The allowed region by the electroweak breaking condition and the 

constraint ml > for k = 1.0. 



In Fig. 2, we find m| < for m 16 < 0.4M and m| < m 2 for 

mi6 < 0.6M. In the universal case with mf 6 = m\ = M 2 /3, the light- 
est superpartner (LSP) is the lightest stau. For other values of k, we obtain 
similar results. 

In the open region of Fig. 2, the mass of the lightest neutral CP-even 
Higgs boson is 90 GeV and the lightest neutralino and the other superpartners 
as well as the other Higgs fields are heavier than 170 GeV. 

Next we estimate the bottom quark mass at the weak scale. The present 
experimental value of the bottom mass includes uncertainties. For example, 
in Ref. [E3fl, it is shown 



m b (M z 



2.67 ±0.50 GeV. 



(16) 



On the other hand, the analysis of the T system [^4j and the lattice result 
[|U give m b (m b ) = 4.13 ± 0.06 GeV and 4.15 ± 0.20 GeV, respectively, f\ 



which translate into 



m b {M z 



2.8 ± 0.2 GeV. 



(17) 



§ See also Ref.Ef 
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For example, the case with k = 0.7 in our model predicts m b (Mz) = 3.4 
GeV at the tree level. However, the large tan f3 scenario, in general, leads to 
large SUSY corrections, i.e, m b = \ b (Hi)(l + A b ). Dominant contributions 



to Ab are given [21 



A b = ^M^tan/3 /(m^mf^Mj) 



3tt 



+ ^A t/ jtan(3 J(m|,m|,/i 2 ) (18) 

where M g , m~ bi and are the gluino, sbottom and stop eigenstate masses, 
respectively. The integral function I(a, b, c) is given by 

ab\n(a/b) +bc\n(b/c) + acinic I a) 

I(a, b, c) = — r . 19 

(a — b)(b — c)(a — c) 

The function I (a, b, c) is of order 1/m 2 ^ where m max is the largest mass in 
the particles running in the corresponding loop. The first term of R.H.S. in 
eq.flrgl) is expected to be sizable. Since the tree level predicted value, m b = 
3.4 GeV, is larger than the values given in (|T6| ) and (|T7l), SUSY corrections 
should be negative. That corresponds to fi < 0. This region is also favorable 
for the constraint due to the b — > decay because this region, fi < 0, can 
lead to smaller branching ratio in the large tan f3 scenario than the prediction 
by the SM P7| . Thus we consider only the case with /i < 0. Fig. 3 shows 
prediction of m b (Mz) including the correction for k = 0.7. The curves in 
the figure correspond to m b (Mz) =2.1 GeV and 2.6 GeV, which are lower 
bounds given in flTBj ) and (|17D, respectively. We have small SUSY corrections 
| A;, | in two regions, where M is larger than and mi6 is much larger 
than M. These regions for fi < lead to the large bottom quark mass, e.g. 
m b {Mz) > 2.6 GeV. This behavior is easy to see since | A&| is suppressed when 
M > nij; or M < mj due to the factor I(m^,m^, M~). Also dotted lines in 
the figure show the boundaries for m| < and m| < m 2 , which are almost 

same as those in Fig. 2. The constraint m| x < excludes the region with 
m b (Mz) > 2.6 GeV for m 16 < M < 1.5 TeV. Furthermore, the stau is the 
LSP in the region with m b (M z ) > 2.6 GeV for 1.5TeV < m 16 < M < 3 TeV. 
To realize m b {Mz) > 2.6 GeV and the neutral LSP, it is needed that M > 3 
TeV or M < mi^ [J The region with M can be more constrained 

more precise prediction of the bottom mass, it is necessary to take into accuont 



8 



by the requirement that the LSP should not overdose the universe, i.e., 
^x^ 2 — 1> m the case where sfermions of first and second families are also 
much heavier than the gauginos |3tJ. A bigger value of k, e.g. k = 1.0 or 
1.4, leads to larger SUSY corrections | A 6 | and predicts smaller values of the 
bottom mass for /i < 0. 
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Fig. 3: Predicted values of m^[Mz) for k = 0.7. 



We have discussed the case that squarks and sleptons in the third family 
have the same soft scalar mass mi6 at Mx, which is required by unbroken 
50(10) gauge symmetry. However, if a gauge symmetry breaks reducing 
its rank like 5*0(10) — > GgMi additional sizable contributions to soft scalar 
masses can appear at the breaking scale, which is called D-term contributions 
[f|, These D-term contributions are, in general, proportional to quantum 
numbers of broken diagonal generators. If we specify the model, one can 
calculate their magnitudes pl[ . Here we study the effect on parameter space 



keeping them free parameters. The soft scalar masses at Mx are written as 



171/ 



rrir 



m 16 - m D , 



rrir 



m 



Hi 



M 2 - 2m 2 16 - 2m 2 D , 



m 



H 2 



= ra 

M 2 - 2m 



TL 

2 o™2 
16 



m\ % + 3m 2 D i 



+ 2m 



D 



(20) 



SUSY threshold corrections j28j. For example, the quasi fixed point of the bottom Yukawa 
coupling as well as the top coupling is raised due to SUSY threshold effects in most of 
cases H. 
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in the presence of the D-term contributions, Qiom 2 D , where Qio denotes a 
broken diagonal charge up to a normalization factor. Note that the above soft 
masses satisfy the sum rule (0) even with taking into account D-term contri- 
butions. Because D-term contributions are proportional to broken charges 
and these charges should conserve in allowed couplings. 

A positive value of m 2 D is unfavorable for successful electroweak breaking 
and the constraint m| x > since it reduces the values of m 2 Hl and Tri~ R . 
Recall that m f „ is dominant to m f , with the initial condition m~ = m~ . 
For example, in the case with m 2 D = M 2 /3, the whole parameter space of 
(M, mio) is ruled out from the electroweak symmetry breaking condition. 
On the other hand, a negative value of m 2 D is favorable to increase both 
m 2 Hl and m 2 R , while it reduces m 2 L . The region excluded by the electroweak 
breaking condition easily becomes narrow for m 2 D < 0. The constraint due 
to the lightest stau is also relaxed as m 2 D decreases from m 2 D = 0. However, 
that reduces m 2 L and below a critical value of m 2 D the value of m 2 L becomes 
dominant to m^. Thus below such a critical value of mjj, the allowed region 
becomes narrow due to the constraint for m% . Around m 2 D /M 2 = —0.1 we 
obtain the widest allowed region. Fig. 4 shows the case with m 2 D /M 2 = —0.1 
and k = 0.7. In this case, the constraint due to the electroweak breaking is 
less important, and actually the excluded region by the electroweak breaking 
is out of the region shown in Fig.4. Dotted lines correspond to boundaries 
for m 2 < and m 2 < m 2 . The region with m 2 < becomes narrow. The 
predicted values of the bottom quark mass is shown as curves corresponding 
to 2.1 GeV and 2.6 GeV in Fig. 4. For M > m 16 , we have large SUSY 
corrections |A&| compared with the m 2 D = case. 
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Fig. 4: The excluded region by the constraint m 2 > and predicted values 
of m b (M z ) for k = 0.7 and m 2 D /M 2 = -0.1. 

If m 2 D /M 2 < —0.1, wider region is excluded by the stau mass constraint. 
On top of that, SUSY corrections |A&| become large in the region M > miQ, 
although the opposite region with M <C m\% leads to slightly small SUSY 
corrections compared with the m^, = case. For example, almost half of the 
parameter space is forbidden in the case with m 2 D /M 2 = —0.3 by the stau 
mass constraint. Thus we can not obtain m^Mz) > 2.6 GeV for M > m^, 
while the region with m&(M^) > 2.6 GeV and M m 16 becomes slightly 
wider compared with the mo = case. Therefore small values of are 
favorable for realistic models. Our result is expected to agree with that in 
SUSY SU(5) GUT models with a large tan/3 because the soft scalar mass 
spectrum in the presence of D-term contribution in SUSY 5O(10) GUT 
models has the same pattern as that in SUSY 577(5) GUT. In fact, our 
result is consistent with that in Ref. |i~l| . 

To summarize, we have considered phenomenological implications of re- 
lations (|) and (D within the framework of SUSY 5O(10) GUT. We have 
investigated constraints due to successful electroweak symmetry breaking and 
the positivity of the stau mass squared. These forbid the parameter region 
with small gaugino mass and m 16 < 0AM. We have also estimated the 
bottom mass in the allowed region. Further, the allowed regions are more 
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constrained by other requirements, e.g., the lightest superparticle should be 
neutral and this particle should not overdose our universe. Also we have 
taken into account D-teim contributions to soft scalar masses. Small D-term 
contributions are favorable to realistic models. It is an interesting subject 
to construct a realistic SUSY SO (10) GUT with sum rules refering to our 
result .0 
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